Manifestations of Broken Symmetry: The Surface Phases of Ca 2 - x Sr x Ru0 4 
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The surface structural phases of Ca2-c C Sr cc Ru04 are investigated using quantitative Low En- 
ergy Electron Diffraction. The broken symmetry at the surface enhances the structural instability 
against the Ru06 rotational distortion while diminishing the instability against the RuC>6 tilt dis- 
tortion occurring within the bulk crystal. As a result, suppressed structural and electronic surface 
phase transition temperatures are observed, including the appearance of an inherent Mott metal- 
to-insulator transition for x = 0.1 and possible modifications of the surface quantum critical point 
near x c ~ 0.5. 

PACS numbers: 68.35.Bs, 61.14.Hg, 74.70.Pq, 71.30.+h 



The discovery of superconductivity in Sr2Ru04 cre- 
ated a flurry of experimental and theoretical activity pQ. 
The structural similarity with La2Cu04, the parent com- 
pound of the superconducting cuprates, combined with 
the nonconventional p-wave superconducting order pa- 
rameter makes S^RuC^ a focus of intense investiga- 
tion [2 j. The substitution of Ca 2+ for Sr 2+ yields a 
phase diagram similar to the high-T c cuprates thus of- 
fering another opportunity to study the ground state 
evolution from an antiferromagnetic Mott insulator to 
a superconductor [3j 2J [5]. One advantage of the 
Ca2- x Sr ;E Ru04 [CSRO] compounds is isovalent substi- 
tution between Ca 2+ and Sr 2+ which alters structural, 
electronic and magnetic properties by tuning lattice dis- 
tortions. Numerous theoretical and experimental works 
reveal the intricate coupling of the Ru06 structural dis- 
tortions with the electronic and magnetic degrees of free- 
dom |3l SI El [71 El [91 HOI HI]. Another advantage is 
that CSRO is a layered perovskite compound, thus its 
crystals are ameanable to cleaving. As such, the study 
of a pristine [0 1] surface is possible through in situ 
cleaving under ultra high vacuum conditions, thus allow- 
ing an opportunity to investigate the intricate coupling 
between structure and other active degrees of freedom in 
an environment of broken symmetry. In this work, the 
surface structural phases are determined by quantitative 
analysis of Low Energy Electron Diffraction (LEED I-V) 
spectra and compared to bulk studies p~2j [13] . 

The bulk structural phases of the CSRO family have 
been previously determined by x-ray and neutron scat- 
tering utilizing both powder and single crystal sam- 
ples [5] El HU [15]. Starting from the highly symmet- 
ric I4/mmm symmetry (no RuOe tilt or rotation) of 
Sr2Ru04 shown in Fig. [lj the smaller Ca 2+ cation shrinks 
the unit cell volume while the Ru06 volume remains 
fairly constant. The shrinking cage surrounding the oc- 
tahedron induces a chemical pressure rotating the Ru06 
into an IA\/acd symmetry while maintaining a uniform 
octahedral shape and volume for 0.5 < x < 1.5. When 




FIG. 1: Left: Bulk U/mmm structure of Sr 2 Ru0 4 . Top Mid- 
dle: LEED pattern for x — 0.5 showing p4gm plane group 
symmetry. Glide lines (on which fraction spots are extinct) 
are emphasized by red lines. Top Right: LEED pattern for 
x = 0.1 showing pg symmetry with only one glide line. Bot- 
tom Right: Structural parameters used in describing bulk and 
surface geometries. 



x < 0.5, an octahedral tilt is induced entering into an 
orthorhombic Pbca symmetry. For 0.2 < x < 0.5, a 
temperature (T)-dependent second order phase transi- 
tion is observed with no hysteresis [5]. For Sr 2 Ru04, 
the system instability against the rotational distortion 
is illustrated by a softening of the RuOe rotational £3 
phonon mode [16]. A similar structural instability for 
x ~ 0.5 is characterized by a softening of the RuOe tilt- 
ing H4 phonon mode [17 . Both the 74 1 /acd and Pbca 
phases can be viewed as arising from the freezing of the 
E3 and £4 modes respectively. For x < 0.2 the system is 
always found in the Pbca phase [5] . Across the metal-to- 
insulator transition (MIT) for x < 0.2 a structural phase 
transition is encountered described by a flattening of the 
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Ru06 and larger lattice distortions. While the bulk sym- 
metry does not change across the MIT, the Ru-0(2) oxy- 
gen bond lengths decrease ~ 0.05 A while the Ru-O(l) 
bond lengths increase ~ 0.05 A. In addition, the tilt of 
the Ru06 increases ~ 5° on average [5]. The structural 
distortions yield a smaller c/a-axis ratio in the insulating 
phase while the volume of the Ru06 increases ~ 3%. 

High quality single crystals were grown using the op- 
tical floating zone technique. All crystals were well char- 
acterized and concentrations verified by energy disper- 
sive x-ray analysis. Crystals were cleaved and measured 
in situ with a base pressure of 8 x 10 -11 torr reveal- 
ing pristine [001] surfaces with large micrometer terraces 
observed by STM. While it has been shown previously 
that the surface of Sr2Ru04 reconstructs to form a lower 
symmetry [18], for x < 1.5 the crystals reveal a p(l x 1) 
surface as shown in Fig. 1. All available beams were 
collected at normal incidence and symmetrically aver- 
aged yielding 16 nonequivalent beams for x = 0.1 and 11 
nonequivalent beams for 0.2 < x < 2.0. Total I — V en- 
ergy ranges varied slightly from surface to surface but all 
I—V sets were > 3000 eV with the majority being > 3700 
eV. Theoretical I — V curves generated for structural re- 
finements employed a modified version of the SATLEED 
program described elsewhere [I9j [20] . Due to the glide 
plane symmetry, simulated annealing optimization algo- 
rithms were written taking advantage of bulk space group 
symmetry generators tailored for each surface [21]. In 
addition, the performance of the simulated annealing al- 
gorithms was checked by manual grid searches for a few 
concentrations. Additional fit parameters were included 
to account for possible asymmetric c-axis displacements 
that do not destroy the observed p(l x 1) LEED pattern. 
The Pendry reliability factor (R p ) was used as a measure 
of agreement between theory and experiment [22 . For 
all surfaces studied the refined surface structures yielded 
0.19 < R p < 0.28 indicating excellent agreement between 
theory and experiment. 

All 0.2 < x < 2.0 samples cleaved at room tempera- 
ture (RT) exhibit a p4gm plane group symmetry. The 
glide lines presented in Fig. [I] are due to the rotation of 
the Ru06 about an axis parallel to the c-axis. While 
the expected symmetry for a bulk terminated IAi/acd 
surface (0.2 < x < 1.5) is p2gg, multiple terrace termi- 
nations generate the p4gm symmetry [20 . For x < 0.2, 
a pg plane group symmetry is revealed, also shown in 
Fig. [T] reflecting the symmetry of the bulk terminated 
Pbca structure. The Pbca symmetry is generated from a 
rotation plus a tilt of the Ru06- The tilt destroys one of 
the glide lines and thus only one is evident in the LEED 
pattern shown in Fig. [I] 

The surface structures for the entire series at RT have 
been determined and the results are presented in Fig. [2] 
For x > 1.5, the bulk symmetry is I4/mmm, however, no 
surface analog to the H/mmm symmetry (P^mm) is ob- 
served for any concentration. The surface stabilizes the 
bulk instability against the RuOe rotational distortion [5], 
freezing in the soft zone-boundary £3 phonon mode cre- 
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FIG. 2: (a) Bulk and Surface RuOe rotation angle versus con- 
centration at T — 300 K. The surface equivalent of a bulk ter- 
minated 14/mmm symmetry never exists as a RuOe rotation 
exists for all x. (b) Bulk and Surface RuOe tilt angles versus 
concentration. The solid lines (filled squares) represent bulk 
powder data while the dashed lines (open squares) represent 
bulk single crystal data [5] [15]. The critical tilt concentra- 
tion is x = 0.2 for both bulk and surface single crystals. Tilt 
angles are slightly enhance for x — 0.1. Inset shows the in- 
ward Ca/Sr displacement for different x. Vertical lines show 
transition from p4gm to pg symmetry. 



ating a single p4gm phase from 0.2 < x < 2.0. The 
surface RuOe tilt angles at RT shown in Fig. [2}d are 
more akin to bulk trends as no RuC>6 tilt is encountered 
for x > 0.2, similar to bulk single crystal data [5j [15]. 
The tilts encountered for the x = 0.1 metallic phase are 
larger than those values encountered in the bulk metal- 
lic phase, but are smaller than the bulk insulating phase. 
The largest surface relaxation observed on the CSRO sur- 
face involves the topmost Ca/Sr ions where a significant 
inward motion is encountered for x < 1.0 as shown in 
the Fig. [2]d inset^ The RT structure for 0.2 < x < 1.0 
shows a large 0.1 A Ca/Sr inward motion but for x = 0.1, 
where a tilt already exists, the inward motion is only 0.06 
A. A simple electrostatic argument would indicate that 
when the surface is formed the topmost Ca/Sr-0(2) layer 
would be forced down [23 , but the insert in Fig. [2]o shows 
that it is not that simple. The surface buckling increases 
and is intimately tied to the stability of the RuOetilt. 
While one might expect the creation of a surface to ac- 
centuate the system instability against the tilt distortion, 
the observed trend discussed below indicates the RuOe 
tilt is stabilized by the creation of a surface. 

The RT LEED pattern for 0.2 < x < 2.0 is shown in 
Fig. [l] The glide lines of the p4gm symmetry is evident 
by the extinguished (=bft, 0) and (0, ±ft) spots where ft is 
an odd integer. To investigate the surface high tempera- 
ture tetragonal-to-low temperature orthorhombic (HTT- 
LTO) phase transition, crystals were cleaved at RT and 
subsequently cooled. As the Pbca bulk phase boundary is 
traversed the tilting RuOe octahedral destroys the glide 
line symmetry resulting in the appearance of the (ft, 0) 
beams. One would expect the low temperature LEED 
pattern to be similar to that of x = 0.1. However, such 
is not the case as both the (ft, 0) and (0, ft) beams are 
evident in the LTO LEED pattern revealing a pm plane 
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group symmetry. 

Using integrated (0,3) and (3,0) beam intensity at 
Ei = 176 eV as an order parameter, the surface HTT- 
LTO phase boundary is determined. As the system is 
cooled, broad diffuse (0,3) and (3,0) beams become evi- 
dent for 0.2 < x < 0.5, indicated in Fig. [3] and [4] by a 
temperature T* . Such diffuse beams are typical of short- 
range correlations similar to those observed in neutron 
data [5j [15]. In contrast to neutron studies, the beam 
intensity is nearly constant for a considerable tempera- 
ture range indicating the system instability against the 
tilt distortion but never achieving the Pbca phase. As the 
phase boundary is traversed, the beam intensity dramat- 
ically increases and the beam size shrinks as long range 
order is established. The behavior of both sets of beams 
is similar across the phase boundary and beam intensity 
is the only difference as shown in Fig. [3J The normalized 
order parameter intensity across the phase boundary for 
x = 0.3 is shown in Fig.[3^ revealing T c ~ 170 K, some 20 
K below the bulk value [24]. While previous bulk stud- 
ies demonstrate the lack of hysteresis indicating a second 
order nature for the bulk phase transition [5JH7], a ~ 10 
K hysteresis is observed on the surface. The doping de- 
pendence for T c has been evaluated for 0.2 < x < 0.5 
(x = 0.2, 0.3, 0.4 and 0.5) and the general trend is similar 
to x = 0.3: the surface T c is suppressed from bulk val- 
ues and a hysteresis is always observed. T* is typically 
larger than the bulk transition temperature. The general 
behavior of this surface phase transition is displayed in 
Fig.g 

For x < 0.2 a doping-dependent metal-to- insulator 
transition exists. Upon cooling, a Mott transition occurs 
between a paramagnetic metal and an antiferromagnetic 
insulator. In the bulk, the transition is coupled to a 
structural phase transition [3, 5, 6 . While the structural 
phases for larger values of x can be described by rota- 
tions and tilts of a rigid Ru06 octahedron, it is found 
that in the insulating state, the octahedron is flattened. 
The flattening octahedron across the MIT is character- 
ized by a sharp decrease of the Ru-0(2) bond lengths 
and increased tilts. On the surface, previous studies have 
shown the surface MIT T c to be ~ 20 K lower than the 
corresponding bulk value for x = 0.1 and it is imper- 
ative to understand the role of surface structure across 
the phase boundary [23]. Figs, ^jp and [3^ reveal striking 
deviations between surface and bulk behavior, the sur- 
face structure across the MIT does not change. While 
the Ru06 tilt increases, it does not increase to those val- 
ues encountered in the insulating bulk. In addition, the 
Ru-O(l) basal plane and Ru-0(2) apical bond lengths, 
as well as all other structural parameters, remain static 
through the phase transition. A 3.3% increase in Ru06 
volume is encountered due to a ~ 4° increase in Ru06 
rotation on the surface. While it has been argued that 
the structural distortions across the bulk Mott MIT are 
responsible for the electron localization [10] [11] , the sur- 
face MIT is not coupled to any structural phase transition 
and is purely electronic in character, i.e. inherent [23] . 
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FIG. 3: (a) HTT-LTO phase transition order parameter for 
x — 0.3 showing first order phase transition character with 
hysteresis. The solid squares are the integrated Beam (0, 3) 
intensity at 176 eV normalized to the Beam (2,2) intensity 
at 115 eV. The solid circles are the normalized Beam (3,0) 
intensity. Vertical line shows transition from p4gm to pm 
symmetry while arrow shows the onset of the tilt instability 
(T*). (b) Surface RuOe tilt angles for x = 0.1 across bulk 
and surface MITs. The four closed symbols represent four 
different crystal surfaces studied. The bulk data (open sym- 
bols) are shown for comparison with lines as guides to the 
eye. The bulk data is from neutron powder experiments with 
a T c ~ 170 K [5] while T c in our bulk single crystals is 154 
K. The surface MIT T c = 130 K 23 . (c) Surface Ru-0 bond 
lengths for x = 0.1 across bulk and surface MITs. Neither the 
RuOe tilts nor the Ru-0 bond lengths show evidence (within 
experimental error) of a structural phase transition across the 
surface MIT. 



Lower HTT-LTO transition temperatures and the lack 
of a structural distortion across the surface MIT suggest 
the tilt is stabilized on the surface. In addition, LDA 
calculations reveal the inward motion of the Ca/Sr plane 
interferes with the tilting of the RuOe across the x = 0.1 
MIT [23 . The general trend suggests the inward motion 
of the topmost Ca/Sr ions plays a significant role in both 
the static tilt across the MIT for x = 0.1 and the sup- 
pressed HTT-LTO phase boundary for 0.2 < x < 0.5. 
The inward motion of the top Ca/Sr ions creates a com- 
pression stress which interferes with the RuOe- Theo- 
retical calculations suggest a similar surface compression 
should exist on other perovskite material surfaces but 
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FIG. 4: Surface phase diagram for Ca2-xSr x Ru04. The 
dashed T* line is the temperature where a tilt instability is 
revealed by weak diffuse reflections and the solid T c line is the 
p4gm - pm structural phase boundary. The dotted lines rep- 
resent bulk structural phase transitions [5l ITol [171124] . There 
is no solid line between the metallic and insulating phases for 
x < 0.2 as no structural phase transition exists across the 
surface MIT. Red arrows indicate different concentrations in- 
vestigated in this study. The unshaded regions below the solid 
and dashed lines for 0.5 < x < 0.8 are extrapolations based 
on observed surface trends. 

experimental evidence has been lacking [25] [26] . The ob- 
served CSRO surface trends would suggest the x c = 0.5 
bulk quantum critical point (QCP) should be shifted to 
lower x on the surface. However, initial results near 
the QCP reveal the surface phases to be more complex. 
While bulk studies reveal the HTT-LTO T c = 155 K for 
x = 0.4 [17], a significant surface suppression of the Ru06 
tilt is encountered as no evidence for the HTT-LTO tran- 
sition is observed down to 80 K. On the contrary, weak 



diffuse superstructure reflections are evident at ~ 80 K 
(T*) for x = 0.5 on the surface and the HTT-LTO phase 
boundary is revealed at T c ~ 40 K. Extrapolation of both 
T c and T* to zero in Fig.[4]shows that the broken symme- 
try at the surface will most likely displace or even destroy 
the x c = 0.5 QCP at the surface. Further investigations 
are required to fully determine the existence and position 
of the QCP on the surface. 

In summary, the surface structural phase diagram of 
Ca2- x Sr :E Ru04 has been determined and is presented in 
Fig. [4] The RT surface structural phases follow bulk 
trends with the exception that no I^/mmm symmetry is 
observed on the surface for x > 1.5. Significant devia- 
tions between surface and bulk behavior are encountered 
across T-dependent structural phase boundaries. While 
the RuOe rotation is revealed for all x, a large inward mo- 
tion of the topmost Ca/Sr ions interferes with the RuOe 
tilt. As a result, lower surface HTT-LTO transition tem- 
peratures are observed for 0.2 < x < 0.5 and the surface 
Mott MIT T c is suppressed for x < 0.2. In addition, fur- 
ther significant surface deviations from bulk behavior is 
noted as a hysteresis is observed across the surface HTT- 
LTO phase boundary and the structural transitions ac- 
companying the Mott MIT in the bulk are simply nonex- 
istent on the surface. Implications of the inward motion 
of the top Ca/Sr ions on the QCP at x c ~ 0.5 are not 
yet clear as an unexpected HTT-LTO phase boundary is 
revealed on the surface at T c ~ 40 K for x = 0.5. 
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